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Recent literature concerning the use of nonlinear transformations to evaluate numerically certain 
improper integrals of the first kind involves the determination of a transformation function g to im- 
prove the approximation. By approximating a given integrand f by an integrable function f\ and then 
determining an associated g function for/i, a nonlinear transformation may be constructed which will 
yield an improved approximation of the improper integral of/. 
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1 . Introduction 

H. L. (/ray and T. A. Atchison 1 introduced a class of nonlinear transformations to assist in 
the numerical evaluation of improper integrals of the first kind. A transformation from this class 
is completely determined by specifying the integrand of the improper integral and a transformation 
function g satisfying certain mild restrictions. The purpose of this paper is to exhibit a scheme 
for the selection of g which will yield a good approximation to the improper integral. 

2. Selecting the Transformation 

The Generalized G-transform of Atchison and Gray is as follows: 
If/ is continuous on [a, o°) and 

F(t)=\ f(x)dx-*S as t^™, (2.1) 

G[F; g; t] = l-R(t) (2 ' 2) 



then 



where geO and g(t) ^ a on [a, oo) , lim g(t) = oo, and 



R( ^7te(7IF<7> <"> 

One of the basic theorems proved (see footnote 1) concerns the determination of a differenti- 
able function g such that G[F; g; t] = S for all t ^ t {) ^ a for some class of f unction sf: 
THEOREM: A necessary and sufficient condition that G[F; g; t] = S for all t ^ t ^ a is that R(t) 
is constant for t ^ to. 

In general, for a given integrand/, the determination of a function gfor which R (t) is constant 
is difficult. However, suppose / is an integrable function which approximates / Requiring the 
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corresponding R(t) to be constant we get 

A (g(t) )g' (t) = aft (t) for t ^ to. (2.4) 

If y=g(t), then (2.4) becomes 

My)dy=c,f 1 (t)dt (2.5) 

and integration yields 

j /i (y)dy=c l ff l (t)dt + C2, (2.6) 

where C\ and c 2 are constants. Thus, if 

H(y) = j ft(y)dy, (2.7) 

then any differentiate function g which satisfies the functional equation 

H(y) = dH(t)+c 2 fort^to. (2.8) 

will result in 

G[Fi;#;t]=Si for t ^ t . (2.9) 

where 

Fi(t)=\ fi(x)dx~>S 1 ast-^co. (2.10) 

J a 



Since /i approximates/, then C[F It ; #; *] will be an approximation of S. 
To illustrate this procedure, consider 

F(0=-7= J c-**^-»erfc (a) as*->oo, (2.11) 

V7T J a 

where < a < oo. For y ^ 1, the integrand of (2.11) may be approximated by 

/ l(y )= re -^ (2.12) 

and f\ possesses an integral 

H(y)=-\e-v\ (2.13) 

One may determine a function y = g(t) which satisfies (2.8) by choosing Ci = e~ fc2 and c 2 = 0. 
Such a function is 

g(t) = Vt* + W. (2.14) 

Applying (2.2) to (2.11) and utilizing (2.14) we get 
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} 



F(t) 



Vt 2 +k 2 



FiVFTk 2 ) 



G[F;g;t] 



\-e k 



Vt 2 + k 2 



(2.15) 



Theorem 1 (see footnote 1) may be applied to (2.15) to show that it converges to erfc (a) more 
rapidly than either F(t) on F( Vt 2 -\- k 2 ). This is illustrated by the following numerical values in 
the case a = l, k=\: 



t 


F(t) 


Error 


G[F;g;t] 


Error 


1.0 

1.5 
2.0 


0.0 
.12340436 
.15262147 


0.15729921 
.03389485 
.0046777 


0.15023133 
.15655151 
.15745495 


0.00706788 
.00074770 
.00015574 



If we apply Aitken's A 2 -process to the F{t) column 2 the resulting approximation is 0.16168462 
which is in error by 0.00438541. 

A general procedure for solving the functional equation (2.8) is not available. Also, for different 
choices of C\ and c 2 , other solutions of the functional equation may be determined. 

(Paper 74B3-329) 
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